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Ultracold three-component atomic Fermi gases in one dimension are expected to exhibit rich 
physics due to the presence of trions and different pairing states. Quantum phase transitions from 
the trion state into a paired phase and a normal Fermi liquid occur at zero temperature. We derive 
the analytical thermodynamics of strongly attractive three-component one-dimensional fermions 
with SU(3) symmetry via the thermodynamic Bethe ansatz method in unequal Zeeman splitting 
fields Hi and H^- We find explicitly that for low temperature the system acts like either a two- 
component or a three-component Tomonaga-Luttinger liquid dependent on the system parameters. 
The phase diagrams for the chemical potential and specific heat are presented for illustrative values 
of the Zeeman splitting. We also demonstrate that crossover between different Tomonaga-Luttinger 
liquid phases exhibit singular behaviour in specific heat and entropy as the temperature tends to zero. 
Beyond Tomonaga-Luttinger liquid physics, we obtain the equation of state which provides a pre- 
cise description of universal thermodynamics and quantum criticality in three-component strongly 
attractive Fermi gases. 
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I. INTRODUCTION 

The ongoing experimental advances in realizing degen- 
erate quantum gases in low dimensions (H-Q offer a new 
and compelling motivation for the further study of quan- 
tum many-body systems via exact schemes such as the 
Bethe Ansatz (BA) and low energy effective field theory 
Q. Reducing the dimensionality in a quantum system 
can have striking consequences. The one-dimensional 
(ID) many-body systems 0, H| possess unique many- 
body correlation effects which are different from their 
higher dimensional counterparts. These include the phe- 
nomena of spin-charge separation, universal thermody- 
namics and quantum criticality. 

A recent scheme for mapping out physical properties 
of homogeneous systems by using the inhomogeneity of 
the trap Q has been successfully applied to experimen- 
tal measurements on the thermodynamics of interact- 
ing fermions with a wide range of tunable interactions 
(lOl Hlj|. Moreover, further experimental advances with 
ultracold atoms allow the exploration of three-component 
Fermi gases in the entire parameter space of trions, 
dimers and free atoms [12Hl4| | . This provides a promis- 
ing opportunity to experimentally explore universal ther- 
modynamics and quantum critical behaviour of strongly 
interacting Fermi gases with high spin symmetries in ID. 
In this context, the thermodynamics of ID attractively 
interac ting fermions [8[ has been receiving growing inter- 
est dMi]. 

For spin-1/2 fermions with attractive interaction there 
are three quantum phases at zero temperature: the fully 
paired phase which is a quasi-condensate with zero polar- 
ization p, the fully polarized (normal) phase with p = 1, 



and the partially polarized (ID FFLO) phase where 
< p < 1 at zero temperature [HI, [ltl H(| ■ This theoret- 
ical prediction of the phase diagram for ID fermions was 
recently confirmed experimentally by R. Hulet's group at 
Rice University 16!]. In addition, it was recently proved 
[l8j that at low temperatures, the physics of the gap- 
less phase belongs to the universality class of a two- 
component Tomonaga-Luttinger liquid (TLL). However, 
from the theoretical point of view, understanding the 
thermodynamics of multi-component Fermi gases with 
higher spin symmetry imposes a number of challenges 



For multi-component interacting Fermi gases, the 
phase diagrams become more complicated in the pres- 
ence of magnetic fields due to the richer number of quan- 
tum phases. In contrast to the two-component Fermi gas, 
[HI O US] three-component ultracold fermions give rise 
to quantum phase transitions from a three-body bound 
state of "trions" into the BCS pairing state and a normal 
Fermi liquid (26l - [35| . The zero temperature phase dia- 
grams of the BA integrable ID three-component Fermi 
gas with SU(3) symmetry have been worked out from 
the dressed energy equations [2{| [35[ . It was found that 
Zeeman splittings can drive transitions between exotic 
phases of trions, bound pairs, a normal Fermi liquid and 
mixtures of these phases, see Fig.[T] It is thus very worth- 
while to map out such zero temperature phase diagrams 
to the inhomogeneity of the trap at finite temperatures. 

In this paper, we investigate the finite tempera- 
ture thermodynamic properties of ID three-component 
fermions with unequal Zeeman splitting by means of the 
exact thermodynamic Bethe ansatz (TBA) solution. We 
prove that at low temperatures the system behaves like 
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either a two-component or a three-component TLL in cer- 
tain regimes. Exact finite temperature phase diagrams 
are demonstrated for illustrative values of the Zeeman 
splitting parameters. Quantum criticality with respect 
to the specific heat and entropy as the temperature tends 
to zero is discussed. The equation of state obtained pro- 
vides an exact description of the thermodynamics and 
quantum critical behaviour of three-component compos- 
ite fermions which can possibly be tested in experiments 
with ultracold atoms. 

This paper is set out as follows. In Section II, we 
present the model and the exact BA solution. We also 
derive the TBA equations for the thermodynamics. In 
Section III, we derive the low temperature thermodynam- 
ics by the Sommerfeld expansion method. The universal 
multi-component TLL phases are identified. In Section 
IV, we present the equation of state in terms of polyloga- 
rithm functions from which the quantum phase diagrams 
can be mapped out. Concluding remarks are given in 
Section V. Detailed working is given in the appendices. 
The derivation of the TBA equations is presented in de- 
tail in Appendix A. In Appendices B and C, the iteration 
method is used to derive relevant results for the TBA and 
the thermodynamics. 



II. THE MODEL AND THE 
THERMODYNAMIC BETHE ANSATZ 
SOLUTION 

We consider a ID system of N fermions of mass m with 
spin independent (^-function potential interaction and are 
constrained to a line of length L with periodic boundary 
conditions. The fermions can occupy three possible hy- 
perfine levels (|1), |2) and |3)) with particle number N 1 , 
N 2 and iV 3 , respectivel y. T he system can be described 
by the Hamiltonian |36ll37| 
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where we have included the Zeeman energy Ez — 
Sj=i-^ le z(MB!-S)- The spin-independent contact inter- 
action 51D applies between fermions with different hy- 
pcrfine states so that the number of fermions in each 
spin state is conserved. The inter-component interaction 
cjid is positive for repulsive interaction and negative for 
attractive interaction. For simplicity, we define the inter- 
action strengths as c = rugm/H 2 and the dimensionless 
parameter 7 = c/n, where n = N/L is the linear den- 
sity, and set H = 2m = 1. Although these conditions 
appear rather restrictive, it is possible to tune scattering 
lengths between atoms in different low sublevels to form 
nearly SU( 3) dege neracy Fermi gases via broad Feshbach 
resonances [L2Wl4| . 

In the above equation, the Zeeman energy levels e z are 
determined by the magnetic moments \i B and the mag- 
netic field B. By convention, particle numbers in each of 



the hyperfine states satisfy the relation A^ 1 > iV 2 > N 3 . 
Thus the particle numbers of unpaired fermions, pairs, 
and trions are respectively given by N\ — N 1 — N 2 , 
N 2 = N 2 — N 3 and N 3 — N 3 for the attractive regime. 

In order to simplify calculations in the study of popu- 
lation imbalance, we rewrite the Zeeman energy as Ez — 
—H\N\ — H2N2 + Ne where the unequally spaced Zee- 
man splitting in three hyperfine levels can be character- 
ized by two independent parameters Hi — I — e z (n B , B) 
and H 2 = 4( M |, B) - e, with e = £- = i 4(/4,B)/3 the 
average Zeeman energy. Pure Zeeman splitting (equally- 
spaced splitting), i.e. Hi — H 2 , leads to a smooth phase 
transition from trions into the normal Fermi liquid. On 
the other hand, unequally-spaced Zeeman splitting can 
lead to quantum phase transitions from trions to the 
fully-paired phase and to a mixture of pairs and single 
atoms, see Fig. Q] 

The Hamiltonian ([T]) exhibits a symmetry of U(l) x 
SU(3), where U(l) and SU(3) describe the charge and 
spin degrees of freedom. This model was solved long 
ago by means of the nested Bethe ansatz [3(1 137[ . The 
energy eigenspectrum is given in terms of the quasimo- 
mcnta {kj} of the N fermions by 
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which satisfy the BA equations [36|, [3 
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Here j = 1,...,N, £ = 1, . . . , Mi, m = 1, . . . , M 2 with 
quantum numbers Mi = N 2 + N 3 and M 2 = iV 3 . The 
parameters {A^,A m } are the rapidities for the internal 
hyperfine spin degree of freedom. 

In the thermodynamic limit, N, L — > 00 with n finite, 
the sets of solutions {kj}, {A^} and {A m } of the BA 
equations (J3j> are of certain forms, as discussed in Ap- 
pendix A. For attractive interaction the quasimomenta 
{kj} can form two-body and three-body charge bound 
states, which give a natural description of composite 
fermions, and can also be real (29l. |37|. However, the 
rapidities {A^} and {A m } can form complex spin-strings 
characterizing the spin wave fluctuations at finite tem- 
peratures. 

In the thermodynamic limit, the grand partition func- 
tion [1, [H| Z = tr(e _H/T ) = e" G/T is given in terms of 
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the Gibbs free energy 

G = E - fj,N + E Z -TS 

= E-fiN- HiNi - H 2 N 2 - TS, (4) 

where the chemical potential fj,, the Zeeman energy E% 
and the entropy S are given in terms of the densities of 
unpaired fermions, charge bound states, trions and spin- 
strings, which are all subject to the BA equations ([3]). 
The equilibrium states are determined by minimizing the 
Gibbs free energy, which gives rise to a set of coupled 
nonlinear integral equations - the TBA equations for the 
dressed energies Sai^ct — 1, 2, 3), which are derived for this 
model in Appendix A, with final result 

£i(fc) = k 2 -n-Hi+ Tai*ln(l + e- £2/T )(fc) 

+ Ta 2 *ln(l + e- £3/T )(£;) 
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+ T( ai + a 3 ) *ln(l + e- £2/T )(k) 
+ T{a 2 + a 4 )*Hl + e- e3 / T )(k). (5) 
Here the quantity 
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and * denotes the convolution, 

(a*b)(x)= / a(x - y)b{y)dy. 
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The spin string parameters £ n := a^/a n and Cn '■= 
/t„ associated with particle and hole densities of string 

length n in A and A parameter spaces satisfy the string 

TBA equations 
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The functions T mn and S mn are as defined in Appendix 
A. 

In the thermodynamic limit, the pressure p is defined in 
terms of the Gibbs energy (0} by p = —(dG/dL), which 
includes three parts, p^ , p^ and p^ , for the pressure of 
unpaired fermions, pairs and trions, respectively, where 
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Here we have set the Boltzmann constant &b = 1. 

The TBA equations ([5]) are expressed in terms of 
the dressed energies £i(fc), e 2 (k) and £a(fc) for unpaired 
fermions, pairs and trions, respectively. The dressed en- 
ergies are seen to depend not only on the chemical poten- 
tial fj, and the external fields Hi and H 2 but also on the 
interactions among themselves as well as the spin fluctu- 
ations characterized by the spin-strings ©. We clearly 
see that spin fluctuations are ferromagnetically coupled 
to the dressed energies for unpaired fermions and pairs. 
There is no such spin fluctuation coupled to the dressed 
energy of the spin neutral trion states. The TBA equa- 
tions play the central role in the investigation of thermo- 
dynamic properties of exactly solvable models at finite 
temperature. They also provides a convenient formalism 
to analyze quantum phase transitions and magnetic ef- 
fects in the presence of external fields at zero temperature 
391. 



III. UNIVERSAL TOMONAGA-LUTTINGER 
LIQUID PHASES 

The TBA equations ([5]) and ([H} involve an infinite 
number of coupled nonlinear integral equations which 
hinders access to the thermodynamics from both the an- 
alytical and numerical points of view. In the strong cou- 
pling regime, the dressed energies e a (k) with a = 1, 2, 3 
marginally depend on each other. The spin string contri- 
butions to thermal fluctuations in the strong coupling 
regime and at low temperatures, i.e. T <C Hi and 
T -C H 2 , are negligible. In this temperature regime, 
the TBA equations can be sorted as 

e a {k) w ak 2 -A {a \ a = 1, 2, 3, (10) 
in terms of the dressed chemical potentials 

A™ = , + Hi-^-^\ 
\c\ 3|c| 

4(2) _ 2u + -c 2 + Ho - — - — » (2) - — » (3) 
- ^ + 2° +Id2 \c\ P \c\ P 9\c\ P ' 

=3^ + 2c 2 - * „W - ±p( 2 ' - 1 P ( 3 ). (11) 
\c\ 3|c| \c\ 

In this case we can directly calculate the pressure through 
(HI), with result 
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in terms of chemical potential temperature T and ex- 
ternal fields Hi and H 2 . Using Sommerfeld expansion, 
we obtain the pressure p^ at low temperatures, 
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The fields H\ and #2 may drive the system into a num- 
ber of different phases. In order to extract the nature of 
the TLL physics from the low temperature thermody- 
namics, we first consider the phase in which trions, pairs 



and unpaired fermions coexist. In this coexisting phase, 
we can apply Sommerfeld expansion under the condition 
that the effective chemical potentials for trions, pairs and 
unpaired fermions are greater than the temperature scale. 
Iteration with the defining relations 
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leads to explicit forms for the pressure 
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The detailed derivation is given in Appendix B. For the 
total number of particles fixed, i.e. n = ni + 2n 2 + 3ri3, 
the free energy can be written as 



where the ground state energy is given by 
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where effective chemical potentials /i^s are given by 
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In order to see universal TLL physics, we calculate the 
leading low temperature corrections to the free energy F. 
Substituting and p^ into (fT51) . after some lengthy 
calculation, we obtain the leading temperature correction 
to the free energy 
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Eq = —Hxri\ - H 2 n 2 - ^c 2 n 2 - 2c 2 n 3 
and the velocities are 
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The particle numbers n\ , n 2 and n 3 of different bound 
states can be obtained approximately by collecting terms 
up to order l/|c| in the expressions for the effective chem- 
ical potentials in (|19p -(|21 [) at zero temperature, 
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with final result 
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This result shows that strongly attractive three- 
component fermions behave like a three-component TLL 
for the coexisting phase of trions, pairs and unpaired 
fermions at low temperatures. Similarly, we can extract 
the finite temperature corrections to the free energy in 
other quantum phases. For example, in the coexisting 
phase of trions and pairs, we have the same universal 
form 
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where the velocities t> 2 and vz have the same expressions 
as that given in (j2~4")l with ri\ — 0. In the above equa- 
tions, the free energy and the thermodynamics are given 
in terms of the chemical potential and the effective Zee- 
man fields Hi and -ff 2 . The chemical potential is con- 
venient for practical purposes in experiments with cold 
atoms, where the chemical potential is replaced by the 
harmonic potential /i — /iq — ^mu) 2 x 2 . The relation be- 
tween fj, and total particle number n can be obtained 
from (fT4|) . 

Although there is no quantum phase transition in ID 
many-body systems at finite temperatures due to thermal 
fluctuations, we shall show that the TLL phases persist 
for non-zero temperatures, as noted in another context 
401. 



IV. THERMODYNAMICS AT LOW 
TEMPERATURES 

For strong attraction (I7I ^> 1) three-atom and two- 
atom charge bound states can be stable under certain 
Zeeman fields. The corresponding binding energies of 



the trions and pairs are given by e t = h 2 c 2 /m and 
Sh = fi?c 2 /Am, respectively. At high temperatures 
T ~ £t,£b, thermal fluctuations can break the charge 
bound states while spin fluctuations cannot be ignored. 
However, such spin fluctuations coupled to the channels 
of unpaired fermions and the spin-1 charge bound pairs 
are suppressed by large fields Hi and Hi at low tem- 
peratures. In this regime, the spin string contributions 
to thermal fluctuations can be asymptotically calculated 
from the TBA equations (0 and ([5]), see Appendix C. 
We have 
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where J\ = 2p\/\c\ and J 2 = p 2 /|c| effective spin-spin 
interactions and 



In{z) = — I < ' 

V Jo 



' cos(n9)d9. 



(33) 



We also see clearly that there is no such effective spin- 
spin interaction for the spin-neutral trion bound state. 
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Using the formula (| 12[) , we can write the pressure p^ 
in terms of the polylogarithm function, i.e. 



^ = -^T^Li 3/2 (-e^ 
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for a = 1, 2, 3, where the polylogarithm function is de- 
fined as 
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We emphasize that the pressure given by (jM)) provides 
the exact equation of state through iteration with (|3l>)) . 
The thermodynamics and critical behaviour can thus be 
worked out in a straightforward manner in terms of a 
special polylogarithm function. 



A. Phase diagram in the /i — H plane 



We hrst consider quantum phases in the \i — H plane at 
low temperatures. Although there is no quantum phase 
transition in ID many-body systems at finite temper- 
atures, the TLL leads to a crossover from relativistic 
dispersion to nonrelativistic dispersion between different 
regimes, which may persist at some non-zero tempera- 
tures [HI, H(| • The zero temperature phase diagrams for 
fixed total number of particles have been explored ear- 
lier [2^, [35j]. The phase diagrams in the /1 — H plane 
from which quantum criticality and the finite tempera- 
ture phase diagrams can be mapped out are investigated 
here. At zero temperature, the \x — H phase diagrams 
can be worked out either from the dressed energy equa- 
tions obtained from the TBA equations in the limit 
T — > 0, or by converting the critical fields in the H — n 
plane, which were found in [29| , into the /it — H plane or 
directly using the equation of state ((541 with T — > 0. 



We first work out the phase diagram for equally-spaced 
splitting (Hi = H2) at T = through analyzing the band 
filling in the dressed energy equations [29(. Here we find 
that the critical field for the phase transition from the 
vacuum into the fully trionic phase is fi c > — |c 2 . The 
critical field for the phase transition from the fully trionic 
phase into the mixture of trions and unpaired fermions 
is determined by the set of equations 
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(37) 



It seems to be very difficult to get a general expression 
for fi c from the condition (|37l) . except for in the strong 
and weak coupling regimes. Nevertheless, we can extract 
the phase boundary by numerical calculation for arbi- 
trary strong interaction. The critical field for the phase 
transition from the vacuum into the fully polarized phase 
is given by /i c > —Hi. The critical field for the phase 
transition from the fully-polarized phase into the mixed 
phase of trions and unpaired fermions is 
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with Qi 
Fig. Ha) 



\f\x + Hi. This phase diagram is shown in 



The phase boundaries for nonlinear Zeeman splitting 
are obtained in a similar fashion. Indeed we find that 
all zero temperature phase diagrams are consistent with 
the \x-H phase diagrams which are directly plotted from 
the equation of state (|34D with the temperature T = 
O.OOle;,, see Fig. [5] For simplicity, we used A, B and C 
to respectively denote the phases of unpaired fermions, 
pairs and trions. The phases A+B, B+C, A+C and 
A+B+C stand for a mixture of corresponding phases. 

The quantum phase segments in an harmonic trap- 
ping potential can clearly be discerned from the phase 
diagrams in Fig. [5] The phase diagram in Fig. [^a) is for 
pure Zeeman splitting (Hi = -H^)- The multi-critical 
point in the phase diagram in Fig. [21a) is located at 
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(^p, — % i ) at T = 0. It may persist for some non-zero 
temperatures due to the existence of TLL phases. In an 
harmonic trapping potential, the mixture of trions and 
unpaired atoms is at the centre of the trap, whereas the 
unpaired fermions are at the outer wings when the ex- 
ternal field H > 4eh/3. However, for H < 4et,/3 almost 
the whole cloud is the trion phase due to a large bind- 
ing energy of trions. The mixture of trions and unpaired 
fermions might lie in a very narrow strip in the trapping 
centre. 

Quantum phase diagrams for unequally-spaced split- 
tings are very intriguing. In the phase diagram Fig. HJd) 
the Zeeman splitting parameters are H2 — 2H±. In 
this case, the pair phase is energetically favoured. From 
the dressed energy equations we can find that the phase 
boundaries intersect at (-jf 1 , — §*) at T = 0. In an 
harmonic trapping potential, when the external field 
H > ^ff 1 , the centre of trap is a mixture of trions 
and pairs whereas the outer wings are occupied by pairs. 
However, for H < -^f-, the mixture of trions and paired 
fermions lie in a narrow strip in the trapping centre. The 
trions occupy the outer wings. 

More subtle quantum phases can be tuned through 
nonlinear Zeeman splitting, see the phase diagrams in 
Fig. [2Kb) and Fig. Hfc), where the phase diagrams for 
the chemical potential are shown for the illustrative field 
values H2 = 1.24ffi and Hi = 1.3-ffi. The mixture of 
trions, pairs and unpaired fermions can occur in a certain 
setting of Zeeman splitting among the three lowest en- 
ergy levels. The intersection points in the phase diagrams 
can be easily determined through the equation of state 
(fM)) with such settings, but it seems to be more difficult 
to analytically determine the phase boundaries. These 
subtle quantum phases can be mapped out through the 
new scheme proposed in Q from experimental data in 
trapped ID Fermi gases. In order to understand the na- 
ture of such quantum phases, we turn to the examination 
of the specific heat in the T — H plane. 



B. Specific heat and entropy 

The thermodynamics of the system ([T]) can be ana- 
lytically calculated through the equation of state ([Mf. 
All thermodynamic properties then follow analytically 
through the general thermodynamic relations. Accord- 
ing to the formula for the specific heat c v — (9 2 p/9T 2 ) , 
the phase diagrams as revealed by c v in the T — H plane 
can be easily explored for fixed total density. Here the 
specific heat c v is a function of T, \i, H± and Hi. Thus 
the full c v phase diagram would be four dimensional. In 
order to observe the signatures of the TLL, we take two- 
dimensional contour plots for the c v phase diagrams for 
some illustrative values of Zeeman splitting associated 
with Fig. H 

For pure Zeeman splitting, the gapless phase is de- 
scribed by a two-component TLL phase under a crossover 
temperature (lines of squares in Fig. [3]) which indicates a 



deviation from the linear temperature-dependent specific 
heat 

irT/1 1\ 
37l \vi v 3 J 

The trions and unpaired fermions can form an asymmet- 
ric two-component TLL of composite fermions and single 
atoms for temperatures below the lines of squares. How- 
ever, the trion phase C and unpaired fermions phase A 
form two different single-component TLLs which lie be- 
low the left and right lines of triangles, respectively. In 
the single-component TLL phase the other states are ex- 
ponentially small and thus the system is strongly corre- 
lated. 

For unequally spaced Zeeman splitting (Hi = 2H\) 
the zero temperature phase diagram in Fig. ^d) may 
persist for finite T as long as the excitations are close 
to the Fermi points of each Fermi sea. From the low 
temperature phase diagram Fig. [4] we see clearly that a 
two-component TLL of trions and pairs remains in the 
regime B + C. The gapless phase is described by a two- 
component TLL phase under a crossover temperature 
delineated by a deviation from the linear temperature- 
dependent specific heat 

Cv^ — l 1" — ) ■ 40 

3ft V"2 v a J 

In this case a TLL of hard-core bosons of composite 
fermions lies below the right line of triangles. 

For unequally spaced Zeeman splitting (H 2 = 1.2Ht) 
the three-component TLL (A + B + C) and two- 
component TLL (A + B) phases may persist within cer- 
tain regimes in the T—H plane, see the lines of squares in 
Fig. [5j Beyond the universal crossover temperatures one 
of the excitations among the states of trions, pairs and 
unpaired fermions exhibits nonrelativistic dispersion. In 
the three-component TLL phase, i.e. where trions, pairs 
and unpaired fermions coexist, the specific heat is given 
by the linear relation 

ttT ( 1 1 1 \ 

3ft \Vi v 2 V3 J 

We see clearly that the equation of state ([Mf provides 
a precise description of the thermodynamics and critical 
behaviour of composite fermions. 

In Fig. [51 we demonstrate that the entropy exhibits a 
peak as the driving parameter chemical potential varies 
across a phase boundary in the fi — H plane, see Fig.[2jc). 
The entropy curves are shown in Fig. [6] for the indicative 
values Hi = 1.2sb, 1.32ef, and Hi = 1.38ef,. In this exam- 
ple, the chemical potential thus varies across the differ- 
ent phase boundaries in Fig. [He) at which the quantum 
phase transitions occur. The entropy peaks in Fig. [B] are 
located in the phases with higher density of states. 
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V. CONCLUSION 

In conclusion, we have studied the thermodynamics 
of ID strongly attractive three-component fermions in 
the presence of nonlinear Zeeman fields via the thermo- 
dynamic Bethe ansatz solution. The pressure and free 
energy have been analytically calculated in terms of the 
chemical potential /u, temperature T and Zeeman fields 
H\ and H2 for a parameter regime T <C £&, £t, Hi, H% 
and 7 ^> 1. Here et and £t are the binding energies 
for a bound pair and a trion, respectively. This physical 
regime covers the presently accessible experimental pa- 
rameter regime @ . The universal thermodynamics of the 
asymmetric two-component and three-component TLLs 
has been identified at low temperatures. Beyond a cer- 
tain crossover temperature, at least one of the underly- 
ing dispersion relations for the composite particles is no 
longer linear and exhibits rich thermal excitations. 

We have derived the equation of state (|34|) from which 
quantum criticality and quantum phase transitions can 
be mapped out. The equation of state provides the nec- 
essary information to describe the quantum regime near 
quantum critical points. The scaling functions and criti- 
cal exponents can be obtained from the equation of state 



following the approach for the two-component model 
(4l| . With regard to the harmonic trapping of three- 
component fermions, quantum criticality can be mapped 
out through the specific heat phase diagram in the T — fi 
plane. For example, for equally-spaced Zeeman splitting 
with Hi = 1.36sb, the critical behaviour of the system 
can be conceived from the specific heat phase diagram in 
the T—fj, plane, see Fig. [7] Our results thus open the way 
for further study of quantum criticality in ID many-body 
systems via their exact Bethe ansatz solution. In this 
case for systems of three-component ultracold fermionic 
atoms. 
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Appendix A: Derivation of the TBA equations 



For the ID three-component fermion system we consider, there are three kinds of states in the system, i.e., unpaired 
fermions, pairs and trions. In the thermodynamic limit and at zero temperature, there are three kinds of quasimomenta 
solutions to the BA equations ([3]). These are real {ki}, with i = l,...,Ni for the unpaired fermions, complex 
roots {ki — Ai ± ^ i I c| }• with £ = 1, . . . , N2 for bound pairs and three-body bound states {k m = X m ± i|c|, A m } with 
to = 1, . . . , iV 3 for trions. 

For finite temperatures, there are also spin strings for spin rapidities A and A, which are characterized by the 
string- hypothesis 

= Ay-i(n+l-20)i, (Al) 
Xf = A; l -i(n+l-2/3)i. (A2) 

where n is the length of the string, j labels the number of strings of length n, and A™ and A™ are the real parts of each 
A and A string. At finite temperatures, there are N[ real quasimomenta kj, N 2 real A,,- and real Xj. The number 
of A( n )-strings is Mi n and the number of the A^™^-strings is Mm- These quantum numbers satisfy the conditions 



Mi = N' % + 2A^ + ^ nMin, (A3) 

71=1 

OO 

M 2 = N' z + ^nM 2n . (A4) 



n=l 

Substituting these three sets of solutions into the BA equations ([3]) gives 

N 2 , A , -I I in N 3 1 \ .1 I 00 Mi„ 



k . L _ -|-| fcj — A; + i|c|/2 k, - \ + i|c| ^ kj A" + inkj/2 
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for unpaired fermions and 

2iA 3 L = tt A, - M + ijcj / kf ] - A, + ijcj/2 fc} ZJ - A, + ijcj/2 \ ff Ajj^rf 3i|c|/2 A 3 - A, + i|c|/2 

,=1,^ - Ai - w U x) - A i - ^i/ 2 *f ) - k j - 'Mi 2 J li Kj " A < ~ 3i|c|/2 Aj " Xi " i|c|/2 

oo Min n . . n .B . . i i 

n=l !=1 ,8=1 -3 1 I C 'I 

for paired fermions. Notice that Eq. (|A6|) has explicit singularities from the terms in the () bracket. To overcome 
this, we write the second term of equation of ([3]) as 



^ A, -k l+ ijcj/2 f kf } Aj ijcj/2 fcf - A, - i|c|/2 \ 
JL A ^ - k > - { W 2 \kf ] A j + i\c\/2 kf A, + MIV 

oo Mi„ n . . n a ., I oo M 2 „ n . \n,a .1 1 ; n 

-nnn f: + nnn ':^; <«> 

n= l m= l^ =1 Aj-- A ro -1|C| n=lm , =la=1 A J A m' + l \ c \l* 
which shows the spin flipping of A"-strings. Substituting (IA7|) back to (|A6I) gives the revised form 



N' N' N' 

a 2iAj L = TT Aj - h + i|c|/2 y4 Aj - A; + ijcj A,- - A; + 3i|c|/2 Aj - X l + i|c|/2 

1 Aj - hi - i|c|/2 11 Aj - A ( - ijcj 11 A, - A, - 3i|c|/2 A, - A, - i\c\/2 



L-- 



J- 1 , 11 A,- - A} 1 - m|c|/2 

n=l 2=1 1 1 " 

for pairs without singularities. Similarly, the equation for trions is 



3 3iAj l = y) Aj - fc, + i|cj -p- A, - A f + i| c |/2 Aj - A ; + 3i|c|/2 
11 Xj -h- ijcj 11 A, - A, - i|c|/2 A, - A, - 3i|c|/2 



H -V- -V' ■ i'- Aj - At + 2i|c| 

ll Aj --A i -i|c|A^-A ; -2i|cr 1 1 



The BA equations for the spin parts are 



^ A™ - fc, - imjcj/2 _ tt tt tt A™ - A" + (m + n + 2 - 2/3)i|c|/2 A™ - A? + (m + n - 2/3)i|c|/2 
11 A™ - h +im\c\/2 ~~ ~ 11 11 11 A™ - A\ l - (m + n + 2 - 2/3)i|c|/2 A™ - AJ 1 - (m + n - 2/3)i|c|/2 

00 M2 " " A™- Ap + (m + n+l-2/3)i|c|/2 
A" 

n=l (=1 ,9=1 3 



X 11 11 11 A m_ A r l _( TO + n + l_ 2 /3)i| C |/2' { 1 



^ A™ - A ; + im|c|/2 ~ " A™ 1 - X? + (m + n + 2 - 2/3)i|c|/2 A™ - X? + (m + n - 2/3)i|c|/2 

11 A™ - A, - im|c|/2 " ~ 11 11 11 A™ 1 - Af - (m + n + 2 - 2/3)i|c|/2 A™ - AJ 1 - (m + ra - 2/3)i|c|/2 

x ft ff TT *r-A? + (™ + "+l-2fli|c|/2 
„=i fi M Af - A? - (m + n + 1 - 20)i|c|/2 ' 1 j 
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Defining the function 6(x) = 2 arctanx and taking the logarithm on both sides of the above equations gives 



N * U A N * , , oo M ln 

ki L - Mj+Ee( ^ ) + I: ^ ) + I:I:e( ^ 

£=1 11 £=1 11 n=l (=1 1 



— ) 



N[ N' 2 N' 3 r 

z=i 11 i=i 11 i=i 

oo M 2 „ . .„ 



A, -Ak L A - A, . A -A, 



)i=i /=i 



N[ \ -h N ' 2 



£=i 



f)(-l—-±\ J. Of A J " A ' \ 



JV' 



+E 



(A12) 



(A13) 



(A14) 



AT - k n 



oo Mi„ 



A™ - A 



n oo M 2rl 



AT - X r ; 



n=l £=1 



n=l £=1 



(A15) 



AT' 



A" 1 - Af 



oo M 2 „ 



XT - X? . 



oo Mi„ 



a™ - Ar 



E*( r W" ) = 2 ^ n) + EE^(^)-EE^(^r^) 



n=l £=1 



n=l 1=1 



(A16) 



Here c' = c/2 with 7j, Jj, ifj, jj™' integers or half-odd- integers depending on the quantum numbers. 



The functions 9 m „ and c mra are defined by 
9 



^m+n + 2# m+ „_ 2 + • • • 2^| TO _ n | +2 + 0| m - n | > for to 7^ n, 
20 2 + 26*4 + --- + 26» 2 „_ 2 + 2 „, form = n, 



9m+n-l + #m+7i-3 + " " " #|m-7i|+3 + 0|j7i-n|+l> f° r m n, 
h+@3-\ I" #2n-3 + #2ti-1, for to = n. 



Finally we define the functions 



N 2 , > N 3 1 \ 00 M ln . .„ 

1=1 |C| 1=1 ^ |C| n=ll=l n|C| 



(A17) 



i ' ( A) = 2AL-5:«(^)-E«(^r)-E 



£=1 



£=1 



£=1 



£?r A ~ A h I ^ A ~\ 



-EE<^>. < A18 > 

71=1 £ = 1 1 1 



N' r 



«w = 3A^-E^(^)-E 



A - A; A — A; 



^3 

E 



(A19) 



and 



^ A m _ J.71 00 M i« 

j m (A") E * 



£=1 



m\c 



,X m -A 



E E 1/1 ■ ) + E E " m "( \ c i\ 1 )' 

n=l £=1 1 1 n=l £=1 1 1 



00 M 2 „ 



o°^ J\^n 



k ^ m ) - E ^) - E E ^(^) + E E ^(^V 1 )- 



(A20) 
(A21) 



£=1 



71 = 1 £=1 



71=1 £ = 1 
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In the thermodynamic limit, we then define 



^ = 2n(p 1 (k)+p h 1 (k)) 



dj'(A) 

dA 
dk 1 (A) 

dk 
djrJjn 

dA n 
dk n (X n ) 

dX n 



(A22) 

2^( P2 (A)+p 2 h (A)), (A23) 
27r(p 3 (A)+pJ(A)), (A24) 
27rK(A")+^(A")), (A25) 
27r(r„(A")+r^(A")). (A26) 



where pi and p^ for i = 1, 2, 3 are particle and hole densities in fc-space, and er„, er^; and r„, t^ - are particle densities 
and hole densities for strings with length n in A-space and A-space. 
Thus we have the integral equations 

= Pi +Pi +a x * p 2 + a 2 * p 3 + J^a„ * a n , (A27) 

n 

— = 92 + P2 + a i * Pi + a 2 * 92 + (ai + a 3 ) * p 3 + V" a n * r n , (A28) 

n 

= P3 + P 3 + a 2 * Pi + (ai + a 3 ) * p 2 + (^2 + a 4 ) * p 3 , (A29) 
for the particle and hole densities and 

a n * pi = <7 n + Cr£ + ^ T mn * CTm - S mn * T m , (A30) 

m m 

&n * P2 = 7"n ~l~ 7"n ~t~ ^ ^ T mn * T m ^ ^ 5 mn * (7 m . (A31) 

m m 

where the functions T mn and S mn are defined as 



7T 

_3_ 

2^ 



Jmn - \ 2a 2 + 2a 4 + --- 



_ 2 H 2a| m _„| +2 + 0| m _„|, form^n, 

+ 2a 2 „_ 2 + a 2 „, for m = n, 



and 



flm+n-l + + " " " &|m-n|+3 + a |m-n|+l, f° r Wl ^ 11, 

ai+ a a -\ + a 2n -3 + a 2n -i, for m = n. 

The energy per unit length can now be written as 

E/L = J k 2 Pl (k)dk + J (2k 2 - C ^)p 2 (k)dk + J (3k 2 -2c 2 )p 3 (k)dk. (A32) 
The total particle number and magnetic numbers are 

N/L = J ( Pl (k) + 2 P2 (k) + 3p 3 {k))dk, (A33) 

Mi/L = J P2 (k)dk + 2 J P3 (k)dk+ X] n / ^(A")dA n , (A34) 

M 2 /L = J P3 (k)dk + J2 n j T n (X n )dX n . (A35) 
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The entropy per unit length is 

S/L = J (( Pl + p\) ln( Pl + pi) - Pl ln Pl - p\ hxp^dk + J ((p 2 + p\) rxi(/9 2 + p 2 l ) - p 2 lnp 2 - p§ hxp^)dk 
+ / {(Pa + Pl) HP3 + Pa) ~ Ps lnp 3 - pl In p%)dk 
+ ^2 ((°™ + CT «) ln ( cr " + a n) - °n lnpiCT„ - o\ \na^)dk 

n J 

+ 53/ + ln ^n + a «) - CT « i 11 / 9 !' 7 " - a 'n ln^ l )dfc, (A36) 

n ^ 

where Inn! w nlnn has been used. 

The Gibbs energy ^ per unit length is 

G/L = J k 2 Pl dk + J (2k 2 - C —)p 2 dk + J (3k 2 ~ 2c 2 )p 3 dk -p J ( Pl + 2p 2 + 3p 3 )dk 

-Hi Upi-2^2 nfJ n + nT n)dk -H 2 I (p 2 + Y ncr n - 2 ^ nT n )dk - TS. (A37) 

n n n n 

Finally, the TBA equations follow by taking the variation of equation (1A37[) and setting it equal to zero, i.e. SG/L = 0. 
In this way 

oo 

In m = (k 2 -p- Hi)/T + ai * ln(l + r^ 1 ) + a 2 * ln(l + t^ 1 ) - ^ a„ * ln(l + (A38) 

n=l 

In 772 = (2A 2 - ic 2 - 2/j - # 2 )/T + a x * ln(l + ijf 1 ) + a 2 * ln(l + Jfe 1 ) + (ai + a 3 ) * ln(l + 773- 1 ) 

- ^a^lna + C 1 ), (A39) 

n=l 

3A 2 — 2c 2 — 3/i 

In 773 = + a 2 * ln(l + rj^ 1 ) + (a x + a 3 ) * ln(l + t^ 1 ) + (a 2 + a 4 ) * ln(l + r]^ 1 ) (A40) 

and 

ln£„ = ra(2ffi - H 2 )/T + a n * ln(l + e -^/ T ) + £ T mn * ln(l + - £ S mn * ln(l + Q 1 ), (A41) 

m m 

In C„ = n(2H 2 -H 1 )/T+ a n * ln(l + e -^/ T ) + £ T ro „ * ln(l + Q 1 ) - ^ S ™« * ln ^ + C*)> ( A42 ) 

m m 

in which we define 77^ = p\j pi (i = 1,2,3), £„ = <7£/V„ and ( n = T%/r n . The TBA equations (|A38|) - (|A42|) can be 
written in the form ([5]) with £j = Tln?7 a (a = 1, 2, 3). 

Appendix B: Details of the Sommerfeld expansion 

In this Appendix, we show that accurate expressions for p and p can be obtained by iteration after Sommerfeld 
expansion. We first use equation (fTTj) to rewrite the pressure (fTB")) in the form 



,w « A (M ( 1 ))3/ 2 d + 2! (Ml U _ _ J^!M 



3/2 

(Bl) 



P 3vr ^ j H sU 3 »J I 3|c| M (3) 9|c| M (3) 3|c| M (3)J ■ (Bdj 
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We can extract explicit analytic results appropriate for the strong coupling regime |c| ^> 1 from these equations 
by iteration and neglecting higher order terms. Making use of equations (1141) and n = rii + 2n 2 + 3n 3 , we obtain 
expressions for n, m, n 2 and 713. After a lengthy calculation, we then obtain fi' 1 ', /i' 2 -* and in terms of ni, n 2 
and 77.3, namely 



(i) 



/' 



(2) 



16 (ii {2) Y /2 
3|c|tt fiW 



T 



16 



1 — —— 



( T 



2 „2 



24 VA( 2 ) 

( /i (l))3/2 / 



8 VAW 
2 \ 12 

) + Id ?r 



4 (M (3) ) 3/2 
]c]tt 

1_ 24 (l(3)) 



7T 2 / T 



8 V^4 (3) 



1 + 



12 V-4 (1) 



3|c|tt M ( 2 ) 



t y 



16 

3|c|tt 



/' 



(2) 



16 {^f 2 
16 / 



8 U (3) ; 



'/i 



(i) 



TT 2 ( T 



(3) ~ __1 7T-2 



24 VA( 3 ) 

( M (D)3/2 / 



|C| 7T \ 

12 U (2) 7 

2\ 



1 



-(— V 

24 ^aw y 



9|c|tt /i( 3 ) 



1 ^vnr 64 ( M (2))3/2 

+ 8 U (1) / / 27|c|tt m (3) 



/,(3) 



(1) 



7T 

1 + tt: 



3| C | TT 
2\ 



1 24 (a«) 



64 



;,(2) 



9|c| 



1 - 



/ T 



2 



24 VA(2) 



12 VA(2) 



(B4) 



(B5) 



(B6) 



Substituting (|B4I) - (|B6|) and (JTSJ) into (jTTJ) and keeping terms to order l/|c|, gives the explicit form for 



A {1) « n 2 7r 2 ( 1 + —7i 2 + 77"3 



„f21 ™2 2 / 1 4 2 16 

2 v |q q 3 I C I 



32 



,4.13) ~ -A-x 1 + — — m + -r-rn 2 + tt«3 
3 V 3|c| 9|c| |c| 



(B7) 
(B8) 
(B9) 



The explicit form for ^jS a ^ without terms can also be obtained as 



,(!) 



,2„2 



T 



6 V n\ir 2 



2 n 3 



3k 



r 



1 + —n 2 + 77713 

c c 



12 V n 2 7r 2 



4 n 3 
27|c| n 2 



7T 

T 



(BIO) 
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,.(3) ~ ^ 2 J l_ 




3|c| n§ 




3|c 



n 3 



(B12) 



Substitute all of these equations back into (|B1[) - (|B3I) and neglecting higher order terms gives the pressure p given in 

(Hz}. 



Appendix C: Dressed energy equations at low temperature 

For the strongly attractive regime, the TBA equations ([5]) can be expanded as 

£l (fc) a fc a_ Ai _ J ff 1 + ^ (2) + 2 ( 3) _ T y, r an( M_ k , )ln{l + r i )dk ^ (C1) 

|c| 3|c| ^7 c 

E2 (fe) « 2fc 2 -2 M -4-^ + nP (1) + n^ (2) + ^ (3) - T E I a n {--k')\n{l + C l )dk', (C2) 
2 |c| |c| 9|c| c 

e 3 (k) « 3fc 2 - 3/x - 2c 2 + Apd) + 8 (2) + 1 (3) (C3) 

c 3 c c 



with 



ln£ n (A) = n{2Hl H2) + ^a n (A) + ^T m „ * ln(l + ln(l + C" 1 ), (C4) 



InCn 



(A) = " (2g2 r gl) + ^a„(A) +E r -* ln ( 1 + C 1 ) - E S ™ * ln(l + Q 1 )- (C5) 



Here J\ = 2p\/\c\, J2 — P2/\c\ and the function a n (k) has the new form 

-M-h&r (C6) 

Compared with equation (JTUJ) , the last terms of (|C1|) and (|C2I) are string terms for spin waves of unpaired fermions 
and pairs, respectively. From ln£ and ln£, we have 

£„(A) « e "( 2 ^- if2 )/ T e 27r - /ia ^ A )/ T e^ T -*^ 1 e-^ s -*^ 1 , (C7) 

C„(A) « e n(2H 2 -fli)/T e 2,J 2a2 (A)/T g E m T^.C; 1 fi - Ero S m „,C. ( C8 ) 

Neglecting higher order correction terms we finally arrive at the dressed energy equations (f3"2"j). The pressure can 
readily be written in terms of the poly logarithm function using (1321) and (fl2)) . 
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FIG. 1: (Color online) Phase transitions from states of trions into (a) normal Fermi liquid, (b) fully-paired states and (c) the 
mixture of pairs and unpaired fermions. The transitions are induced by linear (a) and nonlinear ((b) and (c)) Zeeman splitting. 
Ellipses denote charge bound states. 
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FIG. 2: (Color online) The fi — H phase diagrams at the temperature T = O.OOlet for (a) pure Zeeman splitting and unequally- 
spaced Zeeman splitting (b), (c) and (d). V denotes the vacuum phase, A denotes the unpaired fermion phase, B denotes the 
paired phase and C denotes the trion phase. The phase boundaries are determined by the equation of state (|34fl which are 
consistent with the phase diagrams determined via the dressed energy equations, see the description in the text. 
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FIG. 3: (Color online) The specific heat c v in the T-H\ plane for pure Zeeman splitting with H\ = Hi and total density n fixed. 
An asymmetric two-component TLL remains within a regime below the line of squares between H c \ < H < H C 2. The TLL of 
spin-neutral trion states and the TLL of unpaired fermionic atoms lie below the left and right line of triangles, respectively. 
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FIG. 4: (Color online) The specific heat c„ in the T-H2 plane for H2 = 2H\ with total density n fixed. An asymmetric 
two-component TLL remains within a regime below the line of squares between H c i < H < H C 2- The TLL of spin- neutral 
trions and the TLL of the composite pairs lie below the left and right lines of triangles, respectively. 
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FIG. 5: (Color online) The specific heat c v in the T-Hi plane for Hi — 1.2H\ with total density n fixed. An asymmetric three- 
component TLL remains within a regime below the line of squares between H c i < H < H C 2- An asymmetric two-component 
TLL remains within a regime below the line of pink squares between H C 2 < H < H c s. A TLL of spin- neutral trion states and 
a TLL of unpaired fermionic atoms lie below the left and right lines of triangles, respectively. 
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FIG. 6: (Color online) Entropy per unit length vs chemical potential with T/sb = 0.005 and H2 = 1.3-ffi for different Hi/eb- 
The entropy exhibits peaks in the phases of higher density of states when the chemical potential passes the critical points, see 
text. 
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FIG. 7: (Color online) Specific heat phase diagram in the T — fi plane for equally-spaced Zeeman splitting with Hi — 1.36et. 
The figure shows how the phase diagram extends out from the zero temperature phase diagram. 



